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Abstract 

Lattice gluon propagators are studied using tadpole and Symanzik improved gauge action in 
Landau gauge. The study is performed using anisotropic lattices with asymmetric volumes. The 
Landau gauge dressing function for the gluon propagator measured on the lattice is fitted accord- 
ing to a leading power behavior: Z{q^) ~ (g^)^'^ with an exponent k at small momenta. The 
gluon propagators are also fitted using other models and the results are compared. Our result is 
compatible with a finite gluon propagator at zero momentum in Landau gauge. 
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I. INTRODUCE AND MOTIVATION 



Quantum Chromodynamics (QCD) is believed to be fundamental theory for strong in- 
teractions in Nature. Almost all information about a quantum field theory is encoded in the 
Green's functions and the knowledge of QCD Green's function is important for the under- 
standing of some of the novel properties of QCD like confinement and asymptotic freedom. 
Gluon propagators have been studied by various means in the literature. The ultraviolet (UV) 
behavior of the gluon propagator can be investigated with perturbation theory. In the infra- 
red (IR) region, however, it has to be treated with non-perturbative methods. One can follow 
the continuum approaches: such as truncated Dyson- Schwinger equations(DSEs) exact 

anck type diffusion equation of stochastic 



2, B, y, 0, y, y, q, q, q, q, q. 



The 



renormalization group equations [2] , the Fokker-P 
quantization [sl, or the lattice QCD approaches 
results from these two approaches were also compared. Some agreement were found, how- 
ever, some important issues remain to be clarified. In this paper, we study the lattice gluon 
propagator using tadpole improved lattice actions on anisotropic lattices. While most of the 
previous lattice studies were performed on isotropic lattices, our results can be compared 
with both the previous lattice results and the results using continuum approaches. 

In the study of the IR region of the gluon propagator with lattice QCD, lattice volumes 
have to be large enough since the minimal momenta is inverse proportion to the spatial 
dimensions of the lattice. Investigations with different lattices imply apparent finite volume 
effects 211], and even the extremely asymmetrical box is not safe 22|. A lattice with large 
extensions in all four dimensions will require substantial computational resources. One way 
to proceed is to adopt a larger lattice spacing. Previous studies on coarse lattices suggest 
that the lattice spacing error is still under control if an improved latti_ce_ action is used 
We adopt the tadpole improved gauge action on anisotropic lattices 



16 



17 



23|. 



18|, which 



has the advantage of less lattice spacing error and can be used to generate coarser but larger 
lattice to reach the deeper IR region. We also choose to use anisotropic and asymmetrical 
lattices to further depress spacing errors and to obtain more low momentum modes. 

The anisotropic lattice is a lattice with different cell spacings on different axes. We set the 
equal spatial spacing = as while the temporal spacing is at = -j^ag with the anisotropic 
ratio .^0 = 5. The anisotropic lattice has the advantage of further reducing the spacing 
error while the drawback of further breaking the four-dimensional Euclidean symmetry. 
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Therefore, the renormahzation effect of anisotropic momenta should be taken into account 
which introduces an additional parameter to measure and complicates the determination of 
the physical scale. 



II. LATTICE FORMULATIONS 

A. Gauge Action 

We use the tadpole improved gauge action on anisotropic lattices: 

where is the usual plaquette variable and Ry is the 2x1 spatial Wilson loop on the 
lattice. The parameter Us, which we take to be the forth root of the average spatial plaquette 
value, incorporates the usual tadpole improvement. The parameter designates the bare 
anisotropy. 

B. Gauge Fixing 

The gluon field {^^(x)} associated with a gauge configuration {[/^(x)} is given by 

A,{x + a,e,/2) = :^[U,{x) - Ul{x)] - ^Tr[U,{x) - Ul{x)] (2) 
We fix the gluon field to Landau gauge 

d^^A^ - (3) 
To realize the gauge fix, we first define the gauge transform 

U^{x)^G{x)U^{x)G{x + fiy (4) 

where 

G{x) = e " (5) 
Then we define the gauge fixing functionals 

X,IJ. 



(6) 



We adopt the steepest descents method to minimize [2C 

= ^Y.^r{c,A,{x) + C2A2{x)) (7) 



X 



2V 

where ci and C2 are constants tuned to ehminate the artifact in O(a^) order. The Ai^2 
are defined with the extremum condition of the functionals 

^-^l IttG/^ -a ttG iTTGf^ -A ttG\V 



Aifx) = 



X-pG 

2 cx 5^ Tr [f/^«(x - 2[i)U^{x - fi) - U^{x)U^{x + fi) - h.c] ^ /\,{x) = 



(8) 



Suj^ix) 



C. The Momentum Space Propagator 

The gluon field in the momentum space can be written as 

AM = '-^{P.iQ) - Ul(-l)] - lTr[U,{q) - Uli-q)]} (9) 
where is the discrete momentum in the periodic boundary conditions: 

qt. = —r, n^ = o,i,...,L^-i (10) 

and U^{q) is the momentum space hnk, in the form of 

U,iq) = J2e"''U,{x) (11) 

X 

In the continuum, the momentum space propagator in Landau gauge has the form of 

Ki'l) = '^'^'('^M. - ^^)D{f) (12) 



and therefore the momentum should be corrected as 251] 



^. = fsin^ (13) 
The scalar function D{q^) can be computed on lattice 

^(^^) = _ - i)v ^ ^^^MM-^)). (14a) 

D{^) = ji;r[^^ll^M^M-^)). ^ = o (Mb) 

where A^c = 3, A^^^ = 4 are the dimensions of the gauge group and the space-time, and 
V = Ylfj^ Lfj_ is the lattice volume. 

Another equivalent function is the gluon dressing function: 

Z{f) = fD{f) (15) 



D. The Renormalisation of Anisotropic Momenta 



On a symmetric hyper-cubic lattice, the continuum Euchdean symmetry is broken down 
to hyper-cubic symmetry. When consider the low-momentum modes, one normally uses 
the hyper-cubic lattice momentum squared defined by: (f = (4/a^) sin^(ag^/2) where a 
being the lattice spacing. On an anisotropic lattice, however, the hyper-cubic symmetry is 
further broken down to cubic symmetry. As a result, when consider low- momentum modes, 
one should redefine the anisotropic momentum squared as: 

^' = E^*' + (^^^o)'' ^ = 1,2,3 (16) 

where the renormalisation coefficient C,r is an additional parameter and can be measured 
by fitting physical qualities. To tree level, it is evident that = with being the bare 
anisotropy. However, when quantum fiuctuations are considered, C,r is in general different 
from ^0 and the difference between the two can be substantial for coarse lattices. 



III. LATTICE SIMULATION 



Using the pure gauge action ([T]), de-correlated gauge configurations are generated. Three 
sets of lattices are used in this study and the detailed simulation parameters are listed in 
Table [H For all lattices, one of the spatial dimensions is twice as large as the other two 
which yields more low-momentum modes. The last two sets have almost equal physical sizes 
but different lattice spacings so that the spacing errors and the finite volume effects can be 
estimated. The lattice configurations generated are then gauge-fixed to Landau gauge from 
which the momentum-space gluon propagators are measured. 



A. Fitting the Infrared Exponent k 

Studies on the dressing function in the Landau gauge using Schwinger- Dyson equation 
(DSE) formalism indicates that the gluon dressing function Z[q^) has a power-law behavior 



in the IR region: 27 1 



Z{q^) (X {q'f^ (17) 

The exponent k describes the IR behavior of the gluon propagator. A k value of 0.5 indicates 
that the gluon propagator is finite, while n values above and below 0.5 yield infinite and 



vanishing gluon propagators, respectively. 

We choose to fit Z^ip) measured from our simulations to estimate the value of k, using 
the MPFIT package which can perform non-linear regression in IDL environment. The 
naive errors are used in these fitting procedures since the correlation between the gauge 
configurations is small. Since the power-law model requires Z{q^ = 0) = 0, data point with 
zero momentum is discarded. 

Since the power-law behavior of Z{q^) is expected to be valid only at small momenta, a 
cutoff g^utofj > should be placed on the momentum modes that are measured to set up the 
appropriate fitting range. Ideally, if we have enough low-momentum data points, the result 
of the fitting should not depend on the artificial cutoff parameter g^utofr > 0- -^^^ there 
were not enough low-momentum data points, the final result will have some dependence on 
the cutoff parameter, which turns out to be case of our study. To solve this problem, we 
have tried two methods in the fitting. In the first method, a series of fittings are performed 
with variable cutoff parameter (^cutoff- The value of k are obtained for each of these fitting 
ranges. The final value of n is determined as a extrapolation in the limit ^cutoff 0- The 
results for the value of k are illustrated in Figure [1], [2] and [3] for each set of lattice samples. 
The horizontal axis in these figures are the momentum cutoff parameters ^cutoff- The data 
points are the fitted value of k up to the corresponding cutoff value. The red curves are 
the linear extrapolations of various values of k obtained at different g^utos towards the limit 
q^utos ~^ 0- The extrapolated values at q^^^^s ~^ the final estimates for k using this 
fitting method. 

In the first row of the Figure HI the data points for the dressing function Z{q'^) are 
plotted together with the corresponding power-law behavior fits (the red curve) using the 
extrapolated values of k. The final extrapolated results for k, are also tabulated in the 
first three rows of Table HTl together with the corresponding / d.o.f .. The values for the 
anisotropy parameter = xio are obtained using the second fitting method to be described 
below. The resulting k values at ^cutoff ~^ consistent with a finite propagator at zero 
momentum, i.e. k ^ 0.5 although more definite conclusions will require more accurate 
simulation results. It is also noted that the n values obtained from lattices with similar 
physical sizes (i.e. L12B19 and L16B22) but different lattice spacing are consistent with 
each other within the errors, showing that the value of n is not sensitive to lattice spacing 
errors after tadpole improvement. The k value from the lattice with a larger physical size (i.e. 
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L16B19) is slightly larger from the smaller lattice results by about one standard deviation. 

In the fitting method described above, the value of /t can be estimated quite accurately in 
the IR region. However, since there are not so many low-momentum data points, the value 
of the renormalized anisotropy can not be obtained accurately. In order to accommodate 
more data points at larger momenta into our fitting procedure, we also adopted a second 
fitting method. In this second fitting method, we express the dressing function as the power- 
law behavior times a polynomial in to a certain order: a 

Ziq^) oc (qY^il + a^f + a^f + . . .) (18) 

Using this method, the gluon dressing functions for the three set of lattices are fitted again. 
The three plots in the second row of Fig. H] illustrate the situation of this fitting. The results 
of this fitting procedure are listed in row 4, 5 and 6 of Table [TTl The value of k obtained 
with this method is also consistent with a finite propagator at zero momentum, but with a 
larger compared with the results from the previous method. However, this method yields 
more robust values for the anisotropy C,r- Thus as a final step, we fix the C,r values obtained 
from this second method and perform the fit of k using the first method. These results of k 
are listed in the first three rows of Table [TTl 



B. Fitting Propagator with Other Models 



There are other ways of parameterizing the dressing function Z[q'^) 26|: 



\2k 



g2 + A2' 

z2\2k 



7 fz2 



^pole 



(g2)2« + (A 



2\2k. 



(19a) 
(19b) 



In these parameterizations, another parameter A is introduced. Obviously, in the small 
region they agree with the previous models. Fitting to these two forms are also performed 
and the corresponding results are given in the last six rows of Table [TTl The comparison of 
the fitted curve together with the data points are shown in the last two rows of Figure [H 
The final results from these models are quite unstable for various lattices with large errors 
to the fitted parameters, both n and A. These results seem to be consistent with previous 



studies 



28 



29 



30l |. However, since the fitted parameters have large errors and the values 
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are not quite consistent on various lattices, no definite conclusions should be made from this 
fitting method. 

IV. CONCLUSIONS 

In this paper, the gluon propagator in Landau gauge are investigated using anisotropic 
lattices in the IR region. Improved lattice gauge actions are utilized which reduces the 
lattice artifacts substantially. We also used lattices with one of the spatial direction is twice 
as long as the other two spatial directions. The largest spatial size we used in this study is 
about llfm, allowing us to have more low-momentum modes into the infra-red region. 

The momentum space gluon propagator measured in our lattice simulations are fitted 
using various models. The exponent k, which characterize the power-law behavior of the 
gluon dressing function in the IR region, is found to be consistent with 0.5. This implies 
that the gluon propagator may have a finite value at zero momentum, in agreement with 
the result using other non-lattice methods. This work also confirms that IR region gluon 
propagator can be well investigated by adopting improved gauge action and lattices with 
large volume and coarse anisotropic spacing. 
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TABLE L Lattice parameters. 



Lattice 


Beta 


Physical size 


Configurations 


16 X 16 X 32 X 80 


1.9 


5.65 X 5.65 X 11.30 x 5.65/m^ 


200 


16 X 16 X 32 X 80 


2.2 


4.47 X 4.47 X 8.95 x AAlfm'' 


207 


12 X 12 X 24 X 60 


1.9 


4.24 X 4.24 X 8.48 x 4.24/m'' 


207 



TABLE II: Results for the gluon dressing function Z{q^) from various fitting models. 



Fit Patten 


Parameters 








K 




A 


xVd.o.f. 


Power-law'^ 


L16B19 


fixed at 0.6135 


0, 


.5124 


±0, 


.0044 




0.2054 


J. VJ W id W 


T;1 2R1 q 

±J L ^ I J 1. Cf 


fixed at 6.541 


n 


4976 


± n 


0097 




2320 


Power-law"^ 


L16B22 


fixed at 0.6667 


0, 


.4917 


±0, 


.0129 




0.2473 


Power-law'' 


L16B19 


0.6135 ±0.0316 


0, 


.5117 


±0, 


.0305 




0.0725 


Power-law'' 


L12B19 


0.6541 ±0.0261 


0, 


.4736 


±0, 


.0292 




0.0834 


Power-law'' 


L16B22 


0.6667 ± 0.0240 


0, 


.4461 


±0, 


.0258 




0.0644 


Zcut 


L16B19 


0.5832 ± 0.0362 


0, 


.6745 


±0, 


.0808 


0.5493 ±0.1294 


0.0782 


Zcut 


L12B19 


0.6379 ± 0.0306 


1, 


.0846 


±0, 


.3258 


0.2127 ±0.0868 


0.1471 


Zcut 


L16B22 


0.6564 ± 0.0299 


1, 


.0808 


±0, 


.3261 


0.1164 ± 0.0475 


0.1278 


Zpole 


L16B19 


0.5874 ± 0.0354 


0, 


.6128 


±0, 


.0372 


0.6913 ± 0.0743 


0.0736 


Zpole 


L12B19 


0.6428 ± 0.0293 


0, 


.7630 


±0, 


.0639 


0.4999 ± 0.0335 


0.1234 


Zpole 


L16B22 


0.6583 ± 0.0286 


0, 


.7612 


±0, 


.0634 


0.2763 ±0.0182 


0.1020 



"Extrapolated to zero momenta from a series of fitting. 
''Fitted with EqUH 
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FIG. 1: Extrapolation of the parameter k from a series of fittings with different ^cutoff using the 
power-law model with lattices: F = 16 x 16 x 32 x 80, /3 = 1.9. 




FIG. 2: Extrapolation of the parameter k from a series of fittings with different ^cutoff using the 
power-law model with lattices: F = 12 x 12 x 24 x 60, /? = 1.9. 
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FIG. 4: Fitting the gluon dressing function to various models, (with Eq. \T7\ Eq. [TSl Eq. I19al and 
Eq. I19bl respectively) 
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